Finite temperature phase diagram of a polarised Fermi condensate 
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The two-component Fermi gas is the simplest fermion system displaying superfluidity, and as 
such is relevant to topics ranging from superconductivity to QCD. Ultracold atomic gases provide 
an exceptionally clean realisation of this system, where interatomic interactions and atom spin 
populations are both independently tuneable. Here we show that the finite temperature phase 
diagram contains a region of phase separation between the superfluid and normal states that touches 
the boundary of second-order superfluid transitions at a tricritical point, reminiscent of the phase 
diagram of 3 He- 4 He mixtures. A variation of interaction strength then results in a line of tricritical 
points that terminates at zero temperature on the molecular Bose- Einstein condensate (BEC) side. 
On this basis, we argue that tricritical points are fundamental to understanding experiments on 
polarised atomic Fermi gases. 



Over the past decade, experimental progress in the 
field of cold atomic gases has resulted in unprecedented 
control over pairing phenomena in two-component Fermi 
gases. The ability to vary the effective interaction be- 
tween atoms using magnetically tuned Feshbach reso- 
nances has already permitted the experimental investiga- 
tion of the crossover from a BEC of diatomic molecules 
to the Bardeen-Cooper-Schrieffer (BCS) limit of weakly- 
bound Cooper pairs of fermionic atoms P, [3, H, 0, H, 0] • A 
natural extension of these studies is an exploration of the 
Fermi gas with imbalanced spin populations, especially 
since this system has a far richer phase diagram than the 
equal spin case. As well as exhibiting a quantum phase 
transition between the superfluid and normal states, the 
polarized Fermi gas has been predicted to possess ex- 
otic superfluid phases such as the inhomogeneous Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) state @, where 
the pairing of fermions occurs at finite centre-of-mass 
momentum, and the deformed Fermi surface state Q. 
The exact nature of the superfluid states for the po- 
larised Fermi gas is still the subject of considerable de- 
bate. However, atomic gases provide an ideal testing 
ground for this system, since the particle numbers can 
be varied independently from all other experimental pa- 
rameters, and pioneering experiments have recently been 
performed [n]> [l2|, El • Contrast atomic gases with 
the case of superconductors, where the magnetic field 
used to generate a spin imbalance (via the Zeeman ef- 
fect) also couples to orbital degrees of freedom. 

In this work, we elucidate the finite temperature phase 
diagram of a polarised Fermi gas. While much insight has 
been gained from previous theoretical studies llj, |l5l, EE 
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far a key ingredient of the phase diagram has been over- 
looked: the tricritical point, at which the phase transition 
between superfluid and normal states switches from first 
to second order. By determining the behaviour of the 
tricritical point as a function of interaction strength, we 



can completely characterise the topology of the phase di- 
agram without recourse to an extensive numerical treat- 
ment. Specifically, we shall focus on the uniform, in- 
finite system, and concern ourselves almost exclusively 
with the phase boundary between the normal and ho- 
mogeneous superfluid states. We will, however, discuss 
the ramifications of the inferred phase diagram for the 
trapped system. 



FORMALISM 

Experiments to date exploit wide Feshbach resonances 
and are thus well described by the simplest single-channel 
Hamiltonian, where the two fermion species interact via 
an attractive contact potential 

H - MT«T - Mify. = ( £k ~ ^ c L c k<r 

k <r=t,J. 
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Here, £k = k 2 /2m/ (we set ft = 1 and ks = 1), V is 
the volume, and we define the chemical potential \x and 
'Zeeman' field h such that = /i + h and = fj,— h. At 
present, only pairing between different hyperfine species 
of the same atom has been explored experimentally, so 
we restrict ourselves to a single mass rrif. The interaction 
strength g is expressed in terms of the s-wave scattering 
length a using the prescription: 
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We also derive the Fermi momentum using the average 
density n/2 = (rif + ri;)/2, so that Uf = (S^n) 1 / 3 . 
Throughout our calculations, we will keep n fixed. 

The full phase diagram is parameterised by just a few 
observables: the temperature T = 1/(3, the interaction 
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strength 1/kpa, and the density difference or 'magneti- 
sation' m = n-f — ni- To determine the position of the 
phase boundaries, we must minimise the mean-field free 
energy density 
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with respect to the BCS order parameter A, where £k = 
e k — A* an d -Sk = \/Ck + A 2 - Such a mean-field analysis 
provides a reasonable description of the zero tempera- 
ture phase diagram, but at finite temperature, it neglects 
the contribution of non-condensed pairs to both the den- 
sity n = —dCl/dfi and magnetisation m = —dQ/dh. 
This contribution is necessary to approach the transi- 
tion temperature of an ideal Bose gas in the molecular 
limit, and can be included in the non-condensed phase 
(A = 0) through the Nozieres-Schmitt-Rink (NSR) fluc- 
tuation correction to the energy [32| 
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This gives an estimate of the effect of pair fluctuations 
on the second order phase boundary (but not the first 
order boundary, where A^O). 



PHASE DIAGRAM FOR THE UNIFORM CASE 

Considerable insight can be gained by first examin- 
ing the zero temperature mean-field phase diagram, as 
shown in Fig. [T] The general structure parallels that 
of the two-channel case found in Ref. Since there 

is a gap in the quasiparticle excitation spectrum of 
the unpolarised superfluid, the superfluid ground state 
will remain unchanged for h < mhik-Ek- We see that 
the m — superfluid line in the inset of Fig. [T] corre- 
sponds to an area in the H/ef versus diagram, 
which expands as increases. A key feature of the 

strong coupling side is that for 1/fcpa > 1 the superfluid 
state is able to sustain a finite population of majority 
quasiparticles. This "gapless" [la, |l7( superfluid phase 
is only stable for a < and it thus possesses only one 
Fermi surface. In the extreme BEC limit, this state is 
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FIG. 1: The zero temperature phase diagram within mean- 
field theory for both Zeeman field h/sF and magnetisation 
m/n (inset) versus interaction 1/kFa. There are four different 
phases: the normal (N) state, the phase-separated (PS) state, 
the ordinary superfluid (SF) and the magnetised superfluid 
(SFm). Above the line }i/ef = 2~ 1//3 , the normal state is 
completely polarised (m/n = 1). The red and black lines 
enclosing the PS state are both first-order phase boundaries, 
while the SFm-N transition is second-order, and the SF-SFm 
transition (green line) is at least third-order. The tricritical 
point is represented by orange circles at 1/kFa = 2.368 with 
H/ef = 6.876 or m/n = 1. 



straightforwardly understood as an almost ideal mixture 
of bosonic pairs and fermionic quasiparticles. However, 
as we move towards unitarity, the bosons and fermions 
begin to interact more strongly, leading eventually to a 
first-order phase transition to the normal state. Here, a 
system with fixed m will undergo phase separation into 
normal and superfluid regions if < m < mg, where 
tun.s denotes the magnetisation in the normal and su- 
perfluid phases at h c , the critical field for the first-order 
transition. In the BCS limit (fi = Ef), h c = A/%/2 which 
is less than the quasiparticle gap, so the superfluid state 
is unmagnetised ms = 0, and phase separation occurs for 
arbitrarily low magnetisation, consistent with Ref. 



14] 



For the moment we neglect the FFLO state, but will re- 
turn to this point later. 

A crucial observation is that the line m/n = 1 to the 
right of the region of phase separation can be thought 
of as a continuous zero temperature transition at which 
the condensate is totally depleted. It is thus natural to 
identify the point on m/n — I where phase separation 
starts as a tricritical point. Indeed a Landau expansion 
of the free energy both confirms this and identifies the 
tricritical point at 1/kpa = 2.368. 

With this background, we now turn to the analysis 
of the fate of the tricritical point when temperature is 
finite, beginning with the mean-field description. It is 
well known that there exists a finite temperature tri- 
critical point in the BCS limit l/kpa — » — oo, which 




FIG. 2: Finite temperature phase diagram as a function of 
magnetisation m/n and interaction 1/kpa. The plane at tem- 
perature T — is the phase diagram in Fig. [T] The yel- 
low line represents the locus of tricritical points calculated in 
the mean-field approximation, while the orange tricritical line 
corresponds to mean-field theory plus pair fluctuations. The 
fluctuation correction breaks down in the unitarity regime 
— f < 1/kpa < I, and is thus shown as a dotted line. The 
slice at l/kpa = —I is based on a mean-field calculation and 
it shows the region of phase separation terminating in a tri- 
critical point (yellow circle) at finite temperature, followed by 
a second-order phase transition from the superfluid to normal 
state. Note that the boundary between the FFLO and nor- 
mal states (blue line) defines a small region of FFLO phase 
confined to the BCS side of the crossover, as explained in the 
text. 



is a natural consequence of having a first-order transi- 
tion from the superfluid to normal state at T = and 
a second-order transition at m — 0. First studied by 
Sarma in the context of superconductivity in the pres- 
ence of a magnetic field [33j, the BCS tricritical point 
is located at (T crit /A, h CTit /A) = (0.3188,0.6061) H, 
where A = JtEf exp [— 7r/2|fcFa|] (i.e. at weak coupling 
all energies scale with A). This corresponds to a magneti- 
sation to = 2v(eF)h cr n, where v{£p) — m^y/Tp/v^ln 2 
is the Fermi surface density of states. To investigate how 
the BCS tricritical point is related to the one at zero 
temperature, we must develop a perturbative expansion 
of Eq. ([2]) for small A and general 1/kpa. Doing so, one 
finds (Fig. [2|) that the tricritical point at m/n = 1 is 
connected to that in the BCS limit by a line of tricritical 
points that passes through a maximum somewhere in the 
'unitarity' regime —I < 1/kpa < I. Moreover, for any 



FIG. 3: Finite temperature phase diagram for the two-channel 
model of a narrow Feshbach resonance, where the coupling 
between open and closed channels is weak: 7 = O.f. The 
effective interaction is parameterised by the detuning 5/ef- 
The colour scheme for tricritical lines is the same as in Fig. [2] 

given value of 1/kpa < 2.368, the (T/sf, m/n) phase di- 
agram is highly reminiscent of the 3 He- 4 He system, with 
m/n playing the role of the fraction of 3 He. This is not 
surprising, as the finite to system corresponds in general 
to a mixture of bosonic pairs and fermionic quasiparti- 
cles. Note that even the gapped superfluid can be mag- 
netised at finite temperature due to thermal excitation 
of quasiparticles. Of course, at to = the transition into 
the superfluid state is second order at any point in the 
BCS-BEC crossover. 

It is interesting to examine how the FFLO phase fits in 
with the basic topology of the phase diagram. In the BCS 
limit, we already know that the point where the FFLO- 
normal phase boundary meets the normal-superfluid 
boundary asymptotes to the tricritical point [34} . As- 
suming that the transition from the FFLO state to the 
normal state is second-order (although Ref. [HI found it 
to be weakly first order, this will make a relatively small 
difference), and performing a mean-field analysis, we find 
that the FFLO point of intersection leaves the finite tem- 
perature tricritical point with increasing interaction (see 
Fig. [2]), leading eventually to the extinction of the FFLO 
phase at kpa = —0.35. Note that although this treat- 
ment is somewhat approximate, as we have taken the 
SF-FFLO boundary to be the same as the SF-N bound- 
ary in the absence of FFLO, the point of intersection 
will coincide with that derived from a complete mean- 
field analysis. Moreover, despite all our assumptions, we 
expect the detachment of the point of intersection from 
the tricritical point and the eventual disappearance of 
FFLO to be robust features, since in the BEC regime we 
essentially have a mixture of bosons and fermions. 

The inclusion of the fluctuation contribution Eq. Q 
is crucial for recovering the extreme BEC limit, where 
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it is clear that the (second-order) transition tempera- 
ture asymptotes to TbecK 771 ) = ^bec (1 — m/n) 2 ^ 3 (with 
Tbec ~ 0.218£p), the ideal BEC temperature of a gas of 
bosons of density = (n — m)/2 and mass 2m/. More 
importantly, we find that fluctuations shift the mean-field 
tricritical line to lower temperatures and magnetisations 
on the BEC side, while leaving the tricritical points on 
the BCS side largely unchanged, as expected. However, 
in a broad region around unitarity, we find that the ap- 
proximation underlying Eq. ((3]) generally leads to non- 
monotonic behavior of m(h), with m(h > 0) < for small 
h. We interpret this behaviour as a breakdown of the 
NSR treatment, yielding an unphysical compressibility 
matrix ~d 2 fl/d[j, (7 dfi, a > that is not positive semi-definite. 

To address this problem, we note that the NSR scheme 
is a controlled approximation when we introduce reso- 
nant scattering with a finite width, with the width be- 
ing a small fraction of the Fermi energy 36]. The sim- 
plest such description is provided by the two-channel 
model 37 , [38| . The two-channel description of scatter- 
ing depends upon two parameters: a detuning 5/sf de- 
scribing the distance from the resonance, and a width 7 
of the resonance measured in units of the Fermi energy. 
The one-channel description is recovered in the 7 —> 00 
limit, while the treatment of Gaussian fluctuations is es- 
sentially perturbative in 7, with in Eq. {3]) being 
replaced with — iuj m + 7r~ 1 (q, icj m ), so in this case 
the NSR treatment is expected to be accurate. The re- 
sulting phase diagram is shown in Figure O The zero 
temperature phase diagram coincides with the result of 
Ref. [19j. With fluctuations accounted for, and for suf- 
ficiently small 7, we now find a well-behaved line of tri- 
critical points spanning the crossover region. We expect 
that the true phase boundary at 7 — -> 00 is qualitatively 
similar. 



IMPLICATIONS FOR EXPERIMENT 

We now discuss the consequences of our results for 
trapped gases studied in experiment. Modeling the 
trapped gas by the local density approximation (LDA), 
the spatial dependence of the density induced by the 
trapping potential V(r) is accounted for by a spatially- 
varying chemical potential /i(r) = fi — V(r), with h kept 
constant. In the fi/h-T/h plane, we thus move on a hori- 
zontal line (see Fig.|4]). At sufficiently low temperatures, 
a trapped gas will consist of a superfluid core surrounded 
by the normal state. The transition between normal and 
superfluid states in the trap can be either second or first 
order, depending on whether T/h is above or below the 
tricritical point. Moreover, as long as the temperature 
is non-zero, we can always find a sufficiently small h so 
that T/h lies above the tricritical point. This leads us 
to a key point: if a trapped gas at a given temperature 
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FIG. 4: Phase diagram at 1/kpa — in the {i/h-T /h plane. 
The red and black lines are first- and second-order phase 
boundaries, respectively. The arrows at constant T/h rep- 
resent the trajectories followed when going from the centre to 
the edges of a trapped gas. The two trajectories correspond 
to two different magnetisations of the gas: one greater and 
one less than the tricritical point h CT n . 



and magnetisation has a first-order transition between its 
normal and superfluid phases, then we will always cross 
the tricritical point by decreasing the magnetisation at 
fixed temperature. 

We emphasise that there are qualitative differences be- 
tween first and second order transitions in a trap: the 
former yields a discontinuity in the density and magneti- 
zation at the phase interface, resulting in a form of phase 
separation as seen in recent experiments [HI EH, 0, EH , 
while the latter possesses a density that varies smoothly 
in space. Therefore, the magnetisation and temperature 
at which a tricritical point is crossed should be detectable 
experimentally. In fact, a critical magnetisation for the 
onset of phase separation in a trap has been observed ex- 
perimentally [ll|], and a calculation by Chevy supports 
the idea that this coincides with crossing a tricritical 
point [39]. In addition, the order of the transition will 
have an impact on experiments that use phase separation 
as a signature of superfluidity [12] ■ 

The presence of a first-order transition in the trap can 
be even more pronounced if the density discontinuities re- 
sult in a breakdown of LDA. Experiments on highly elon- 
gated traps already provide evidence for such a break- 
down 11], and one requires the addition of surface en- 
ergy terms at the phase interface to successfully model 
the trapped density profiles [40(. 

An outstanding issue is the experimental detection of 
the gapless SFm phase. While optically probing the 
momentum distribution of the minority species is one 
promising method for detecting SFm [41| , another possi- 
bility is to study density correlations usingjfor example, 
shot noise experiments as suggested in Ref. 42] . A simple 
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mean- field calculation gives (for the uniform system): 
C u (k 1; k 2 ) = (^(kx^Cka)) - (h^ih^)) 

- S kl ,- k2 ^[l - f(E kl + h)- f(E kl - h)f 

where f(E) is the Fermi-Dirac distribution. At T = 0, 
the result is a 'hole' in the correlation function for mo- 
menta less than the Fermi wavevector of the majority 
quasiparticles. Such a measurement would therefore con- 
stitute both a confirmation of the SFm phase and a vivid 
demonstration of the blocking effect of quasiparticles on 
(+k, — k) pairing. 

In conclusion, we have determined the structure of the 
finite temperature phase diagram of the two component 
Fermi gas, as a function of both interaction strength and 
population imbalance, finding a region of phase separa- 
tion terminating in a tricritical point for general coupling 
in the BCS-BEC crossover. A secondary result of our 
work is the demonstration that the NSR scheme yields 
unphysical results in a broad region around unitarity. 
This is significant, as it is widely viewed as offering a 
smooth, albeit uncontrolled approximation throughout 
the crossover. We emphasize that there is no a priori 
reason to believe in the accuracy of the NSR scheme with- 
out introducing an additional parameter, as we have done 
here. The Ginzburg criterion governing the smallness of 
fluctuation corrections is satisfied in both the BCS limit 
where it takes the form (T c /ef) 2 <C 1, and in the BEC 
limit where fcj?a <C 1 is the relevant criterion. But at 
unitarity the shift in the transition temperature relative 
to the mean field value will be of order sf- At the same 
time the upper critical dimension at the tricritical point 
is three, so we may expect that our results there will be 
little changed. 

Finally, we have argued that these tricritical points 
play an important role in experiments on trapped Fermi 
gases (see, also, the subsequent related work on trapped 
gases at unitarity by Gubbels et al. (iH). Indeed, a recent 
comprehensive study of the temperature dependence of 
the phase-separated state at unitarity has yielded experi- 
mental results consistent with the phase diagram outlined 
here Q. 

We are grateful to P. B. Littlewood for stimulating 
discussions, and J. Keeling for help with the numerics. 
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